Single parameter scaling in one-dimensional localization revisited 
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The variance of the Lyapunov exponent is calculated exactly in the one-dimensional Anderson 
model with random site energies distributed according to the Cauchy distribution. We find a new 
significant scaling parameter in the system, and derive an exact analytical criterion for single param- 
eter scaling which differs from the commonly used condition of phase randomization. The results 
obtained are applied to the Kronig-Penney model with the potential in the form of periodically 
positioned 5-functions with random strength. 
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Introduction The hypothesis of single-parameter scal- 
ing (SPS) in the context of transport properties of disor- 
dered conductors was developed in Ref. ID. It was sug- 
gested that there exists a single parameter, conductance 
g, which determines a scaling trajectory g(L), where L 
is a size of a sample. It was soon understood [0 that 
scaling in the theory of localization must be interpreted 
in terms of the entire distribution function of the con- 
ductance. In order to take the fluctuations of the con- 
ductance into account, it is convenient to consider a pa- 
rameter j(L) = (1/2L) In (1 + l/g) instead of g itself |§. 
The main properties of this parameter are that its limit 
7 = "f(L — ► oo) is non-random || and that it has a nor- 
mal limiting distribution for L 3> 7 _1 @]. SPS in this 
situation implies that the distribution function of g is 
fully determined by the mean value of ^(L). This mean 
value, which is the scaling parameter, is close to the lim- 
iting value 7, provided that L 3> 7 _1 - The parameter 7 
is called the Lyapunov exponent (LE), and its inverse is 
the localization length, li oc — 7 _1 , of a particle's wave 
function. Below the random function j(L) will be re- 
ferred to as the finite size LE. The localization length 
of the Anderson model (AM) have been calculated by 
many authors, see Refs. ||,|]Jf|. SPS presumes that the 
variance, a 2 , of j(L) is not an independent parameter 
but is related to LE in a universal way 



= 7 /L. 



(1) 



This relationship was first derived in Ref. B within the 
so called random phase hypothesis, which assumes that 
there exists a microscopic length scale over which phases 
of complex transmission and reflection coefficients be- 
come completely randomized. Under similar assump- 
tions, Eq. (g) was rederived later by many authors for 
a number of different models. According to the random 
phase model, if the phase randomization length, L^, is 
smaller than the localization length, li oc — 7 -1 , then 
Eq. (|l|) is valid. An analytical derivation of the distribu- 
tion function of LE for the continuous Shrodinger equa- 
tion H and numerical simulations of tight-binding AM 



M showed that the inequality l p h.j ^> 1 holds as long as 
disorder remains weak. Accordingly, SPS was believed to 
violate only at strong disorder, when li oc becomes micro- 
scopic. However, recent numerical simulations for a dis- 
ordered Kronig-Penney model (KPM) |8| demonstrated 
that, contrarily to the existing picture, even when disor- 
der is weak and li oc is large, a strong violation of SPS is 
still possible. 

In this paper we demonstrate, that the condition for 
the validity of SPS based upon the phase randomization 
concept is not accurate. We study an emergence and vi- 
olation of SPS in one-dimensional systems without the 
assumption of phase randomization. We calculate ex- 
actly the variance of the finite size LE for AM with the 
Cauchy distribution of site energies (the Lloyd model) 
and derive Eq. (Q) without any ad hoc hypothesis. This 
calculation also produces an analytic criterion for SPS in 
the form 7 _1 3> l s - The new length scale l s , which is 
different from the phase randomization length, is a new 
significant length scale. We show that even in the limit 
of weak disorder, when the localization length is macro- 
scopic, the states at the tails of the spectrum never obey 
SPS. 

One-dimensional models with off-diagonal disorder 
(random hopping models) represent a special case. These 
models demonstrate a delocalization transition in the 
vicinity of the zero energy state |^|, which results in a 
violation of SPS as well as in other unusual phenom- 
ena. In this paper we focus upon regular one-dimensional 
situations, which include models with diagonal disorder 
as well as random hopping models far away from the crit- 
ical point, E = 0. 

SPS also takes place in the regime of weak localiza- 
tion, which occurs at weak disorder in two and three 
dimensional electron systems, as well as in quasi-one- 
dimensional wires in the limit of a large number of con- 
ducting channels Jul. In this case, SPS manifests itself 
in the form of universal conductance fluctuations jl2],|l3| ■ 
In the weak localization limit the mean conductance co- 
incides with (jL)^ 1 and remains the only significant pa- 
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rameter, which also determines non-Gaussian corrections 
to the high moments of the distribution function Jl4) . 
SPS is violated only for very higher moments, when the 
far tails of the distribution function become important 
jl4],[l5|]. This situation is also out of the scope of the 
present paper, since the weak localization limit can not 
be realized in one-dimensional models. 

Calculation of the variance of the Lyapunov exponent 
for an exactly solvable model In this paper we discuss a 
quantum particle on a chain of sites, which is described 
by the following equations 



"0n+l + tpn-1 - Un^n = 0, 



(2) 



where ip n represents the wave function of the system at 
nth site. The meaning of U n depends upon the inter- 
pretation of the model (|^). AM corresponds to Eq. (j^) 
with U n = —E + e„, where E is the energy of the par- 
ticle, and e n is a random site energy. Another realiza- 
tion of Eq. (Q) is a KPM with the potential given by 
the sum, ~Y^V n 8{x — na), of periodically positioned 6- 
functions with random strengths, V n . In this case U n is 
defined as U n = 2cos(ka) + -^p-sin(fca), where k = \f~E~ 
is an energy variable and a is the period of the structure 
HjuSJ. We assume that the parameters e„ or V n are dis- 
tributed with the Cauchy probability density (the Lloyd 
model) : 



Pc{x) = - 



ix (x- £ )O - £5)' 



(3) 



where to = xq + iT, and Q = xq — iT. Parameters Xq 
and r represent the mean value and the width of the dis- 
tribution, respectively. For AM xq — 0, whereas for the 
KPM xq = Vq. The distribution of parameters U n also 
has the form of Eq. (||), where we now denote Uq as the 
center of the distribution, and Tjj as its width. Specific 
expressions for Uq and Tjj in both AM and KPM can be 
obtained straightforwardly. 

Below we calculate the mean value, 7, and the vari- 
ance, a 2 , of the finite size LE, "f(L), in the limit L ^S> 7 _1 
(further, all lengths are normalized by the lattice con- 
stant, a). Following Ref. || we present average LE as 

7 = (In \z n \) , z n = V>„/-0n-i 

and averaging (. . .) is carried out over the stationary dis- 
tribution of the random variables z n . It turns out that 
z n are distributed according to Eq. (|^) with the complex 
parameter £0 replaced by tst , which satisfies the equation 



U + iT t 



It is convenient to parameterize tst as 
£at =pexp(i(p). 



(4) 



(5) 



In this notation 7 = In \p\. Both p — |£ st | and the phase, 
<p(E) — arg(£ st ), depend on the energy, E. 



The variance o~ can be expressed in terms of 

N-lN—n 1 

ct2 = 72 E E (lnN„|ln|z„ +fc |) + -(ln 2 |z n |)-7 2 . 



n=0 k=l 



(6) 

Therefore, a 2 can be expressed through the two point dis- 
tribution, Pi(z n , z k ), of the parameters z n . A joint distri- 
bution of multiple random variables can be expressed in 
terms of the product of marginal and conditional distri- 
butions. The latter probability distribution of z^, under 
the condition that z n is fixed, can be shown to satisfy the 
following recurrent relation [k > n) 



P(z 



■ n \zk+i) = J P (z n \z k ) P c (z k+1 + z^ 1 ) dz k , (7) 

where Pc is the Cauchy distribution, Eq. (^). The ad- 
vantage of the Cauchy distribution is that the recurrence 
relation (0) can be solved exactly. The solution again has 
the form of the Cauchy distribution, Eq. (^): 

P(z n \z k ) = ^- 1 — -, (8) 



7T [Zk ~ tk-n) (Zk ~ tk- n ) ' 

with complex parameters tk obeying the equation 



fe-i 



U a + iT v . 



(9) 



and T k — Im(£k). The same Eq. (^) determines parame- 
ter tst of the one-point distribution of quantities z n || . 
In the latter case however, one looks for the stationary 
solution of the equation, while the conditional distribu- 
tion requires solution of Eq. (^) with the initial condition 
£0 = z n . With the use of the evaluated two-point proba- 
bility distribution, the correlator (In \z n \ In |z n _|_/-|) can be 
represented as 



(ln|2„|ln|z n+ fc|) = 



Im£ s 



m|z|ln|q(z)| 



dz. 



One can substitute this correlator into Eq. (g) and sum 
it over k without further assumptions. After calculating 
the average (in 2 \z n \) the variance a 2 can be expressed 
through the parameters p and ip. Eq. (|J): 



-7 In 



p 4 - 2p 2 cos(2ip) + 1 
(P 2 ~ I)' 



(10) 



dx arctan 







( COS ip — COS X 



O (l/L 2 



Here we are interested in the limit 7 <C 1, i.e., the local- 
ization length is large li oc S> a = 1. In this limit the pa- 
rameters /3(7, if) and £(7, ip) are equal to f3 ~ 27 sin ip and 
C — 1 + 27 2 , respectively. We also assume that (p < ir/2, 
the case ip > it/ 2 can be handled by the replacement 
ip — > it — ip. Eq. ( [To| ) is the main result of our calcu- 
lations. It presents the asymptotically exact (L — * 00) 
expression for the variance of LE in the Lloyd model. 
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Discussion The behavior of the integral in Eq. (|1 
is governed by the parameter 



k = (<yl s )~ 



(11) 



where l s is a new length scale in the system: 

l s = l/sm<p= |&t|/Im£ st , (12) 
and £ s t is given by Eq. (Q). In the limit l s <C 7, i.e., 

k»1, (13) 



Eq. (O) reduces to 



27/^, 



(14) 



implying validity of SPS. Thus, Eq. (13) represents a true 
criterion for SPS. 



It should be noted, however, that our result, Eq. (14), 
differs from Eq. ([!]) by the factor of 2. This difference 
reflects the well known peculiarity of the Cauchy dis- 
tribution - all of its moments, except for the first one, 
diverge. As a result, neither the approximation of the 
Gaussian white noise ||], nor a weak disorder expansion 
[JtJ used to derive Eq. (Q) can be applied to the Lloyd 
model. We have calculated the variance of LE for AM 
with the Cauchy distribution of site energies using the 
random phase hypothesis, and following the approach of 
Ref. . We found that though a proportionality between 
the variance and LE persists, the numerical coefficient 
differs from both Eq. (@) and Eq. (|l|). This result im- 
plies that the phase randomization model is not valid at 
all for the Lloyd model. More important, however, is the 
fact that SPS holds even when the phase randomization 
hypothesis fails. This is an additional confirmation of 
the fact that the real criterion for SPS is given by the 
inequality (|l3| ) rather then by phase randomization. 

Even though the localization properties of the Lloyd 
model are the same as those of generic models |l7j , one 
might question the generality of the new scale l s and cri- 
terion (|l^) in light of the peculiarity of the Lloyd model. 
In order to confirm a generic nature of applicability of 
Eq. (13), we carried out numerical simulations of KPM 
with rectangular barriers of random widths. Statistics 
and the shape of potential in these calculations are con- 
siderably different from the Lloyd model. The results of 
the simulations are shown in Fig. 1 along with t(k) ob- 
tained from our analytical Eq. ([To]). One can see from 
Fig. 1 that the crossover between different asymptotes 
occurs in the same region for both models. This allows 
us to conclude that the crossover length l s and criterion 
( |l3| ) retain their significance beyond the Lloyd model. 

According to Thouless [^8), the phase ip(E) is propor- 
tional to the integrated density of states (p(E) = irG(E) 
(in the case of KPM ip(E) must be reduced to the interval 




1 2 3 4 5 6 7 

FIG. 1. The function r^^obtained from the analytical 
solution of the Lloyd model Eq. ( |Io[ ) (solid line) and the 
numerical simulations of KPM with rectangular barriers of 
random widths (open circles). Note that though the two dif- 
ferent models give slightly different asymptotes of the t(k) 
function, the crossover length is the same for both models. 

[0, it}). For the states close to the center of the initial con- 
ductivity bands, <p(E) ~ ir/2 and consequently l s ~ 1. 
Inequality ( |l3| ) reduces then to Z ioc » 1, which essen- 
tially implies that disorder is weak. However, as soon 
as the energy approaches a band edge, / s grows signif- 
icantly, so that k 1 can coexist with weak disorder, 
hoc S> 1, when (^<C 1 or 7r — ip 1. These inequali- 
ties hold for both AM and KPM at energies outside the 
initial spectrum, i.e., at the tails of the density of states 
for AM and within former band-gaps for KPM. In this 
case l s {E) can be expressed in terms of the number of 
states between E and the closest fluctuation spectrum 
boundary. For AM with the Cauchy distribution these 
boundaries are at ±00. For KPM their role is played by 
energies at which (p(E)/Tr is an integer. The states in 
these regions arise due to rare realizations of the disor- 
der, and can be associated with spatially localized and 
well-separated structural defects. The length l s can then 
be interpreted as an average distance between such de- 
fects. In view of this interpretation of the transition 
between two types of scaling regimes occurs when the av- 
erage distance between these defects becomes comparable 
with the localization lengths of the respective states. 

The variance a 2 in the limit opposite to Eq. ( |T3] ) 
can be conveniently presented in terms of the parame- 
ter t = a 2 L/(2~f): 



(15) 



It is, thus, determined by the scale l s rather than the 
localization length li oc . This equation describes the tran- 
sition from SPS to scaling with two independent param- 
eters hoc and l s . The particular form of the function 
t(k) is probably model dependent. At the same time our 
results suggest that n is a universal parameter, which 
naturally describes the crossover between different types 
of statistics in the localization problem. 

The previous criterion for SPS based upon the 
phase randomization length effectively restricts the 
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strength of disorder The new criterion, Eq. (13), 

separates all the states of the system with a given 
strength of the disorder into two groups, which demon- 
strate distinct scaling properties. The boundary between 
the groups, defined by the condition k = 1, coincides with 
the boundaries of the initial spectrum \Uo\ = 2, provided 
that Tjj <C 1. For the states inside the original spectra, 
SPS holds as long as the disorder remains weak. The 
length scale l s in this case is of the order of the lattice 
constant a and, thus, has no scaling significance. Nev- 
ertheless, it is important to notice that this length scale 
differs from the phase randomization length. Indeed, it 
was found numerically in Ref. JtJ that l p h increases to- 
ward the center of the band E = as 1/E. Absence 
of the phase randomization for E = was also obtained 
analytically in Ref. p9[. Both numerical, Ref. Jjj, and 
analytical, Ref. calculations suggest that SPS still 
holds at E = 0, in spite of the fact that the phase does 
not randomize and the standard weak-disorder expan- 
sion for LE fails (|||at this energy. At the same time, 
applying criterion ( |13[ ) one finds that SPS in this case 
violates only when disorder is so strong that the localiza- 
tion length becomes of a microscopic magnitude, li oc ~ 1 
m accord with Ref. §. Indeed, l s ~ 2/^4- J7 2 and de- 
creases toward the center of the band, where it assumes 
the minimum value l s = 1. This allows us to conclude 
that initially conducting states always demonstrate SPS 
in a meaningful scaling regime l\ oc 1. 

The situation for the states from the former band-gaps 
is totally different: SPS violates when both lengths, l s 
and /; oc , are macroscopic. Such a violation is significant 
from the scaling point of view. In the case of KPM Z/ oc 
remains macroscopic throughout entire band-gaps, pro- 
vided that the energy is sufficiently high, k ^> Vq. Con- 
trarily, h oc ~ l for those states of AM, that are not very 
far from the boundary. A strong violation of SPS in this 
case coincides with the system being driven out of the 
scaling regime. 

Change in the strength of disorder affects the two scal- 
ing lengths li oc and l s differently. In the case of the 
gap states, li oc only weakly depends upon disorder. It 
is approximately equal to the penetration length, which 
would describe tunneling through the system in the ab- 
sence of disorder. Disorder related corrections to this 
quantity are of the order of Ty. The parameter Ij , on 
the contrary, is linearly proportional to Tjj. Increase in 
disorder, therefore, causes the critical parameter k to in- 
crease. Thus, the system can crossover between k <C 1 
and SPS (k ^> 1) regimes with an enhancement of the 
disorder. The results seem paradoxical since the restora- 
tion occurs because of increasing disorder, which must, 
however, remain small enough to keep the system within 
the scaling regime. This effect is also more important for 
KPM than for AM, and was observed numerically in Ref. 
H for a periodic-on-average random system. 

Another interesting feature, which is specific for KPM 



and absent in AM, is a presence of resonance states where 
7 = regardless of the disorder. Though this feature 
is unstable with regard to a violation of the strict pe- 
riodicity in the position of the 5— functions, it is also 
present in some other models, such as random supcrlat- 
tices, and deserves a discussion. At the resonance, both 
7 and I J 1 vanish. Their ratio, k, however, remains finite 
and experiences a jump at the resonance point from the 
value of 2Vo/r 3> 1 at the band side of the resonance to 
r/2Vo <C 1 at the gap side. It means, that the change of 
statistics of LE at the resonance states from SPS behav- 
ior at the band side to two parameter scaling at the gap 
side occurs discontinuously. 

We are indebted to A. Mirlin for a useful discussion. 
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menting on the manuscript. Work at Queens College 
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